Robust Ising Gates for Practical Quantum Computation 
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I describe the use of techniques based on composite rotations to combat systematic errors in 
controlled phase gates, which form the basis of two qubit quantum logic gates. Although developed 
and described within the context of Nuclear Magnetic Resonanace (NMR) quantum computing 
these sequences should be applicable to any implementation of quantum computation based on 
Ising couplings. In combination with existing single qubit gates this provides a universal set of 
robust quantum logic gates. 

PACS numbers: 03.67.-a 



Quantum computers JIJ use quantum mechanical ef- 
fects to implement algorithms which are not accessible to 
classical computers, and thus are able to tackle otherwise 
intractable problems They are extremely vulnerable 
to the effects of errors, and considerable effort has been 
expended on correcting random errors arising from deco- 
herence processes ||, ||, ||. It is, however, also important 
to consider the effects of systematic errors, which arise 
from reproducible imperfections in the experimental ap- 
paratus used to implement quantum computations. 

It makes sense to consider systematic errors as some of 
them can be tackled relatively easily. In the context of 
Nuclear Magnetic Resonance (NMR) quantum computa- 
tion H [t], ||, [)[ |uj, |ll| the concept of composite rotations 
(also called composite pulses) [l^, |l4| , |i~5| has been 
used to tackle both off-resonance effects ]16[ and more 
recently pulse length errors JT7[. All these results, how- 
ever, have been concerned with single qubit gates, and 
while these play a central role in quantum logic circuits 
it is also vital to consider two qubit controlled logic gates, 
such as the controlled-NOT gate. 

The combination of a complete set of single qubit gates 
and the controlled-NOT gate is universal for quantum 
computation Jl^j , that is any quantum logic circuit can 
be built entirely from single qubit gates and controlled- 
NOT gates. There are, however, other two qubit gates 
which can be used to form a universal set, most notably 
the controlled phase-shift (controlled-er 2 ) gate [|To). This 
performs the transformation 
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while leaving other eigenstates unchanged, and so is de- 
scribed by the matrix 



(2) 



Applying Hadamard gates to the target spin before and 
after the controllcd-cr z gate converts it to a controlled-o^ 
gate, that is controlled-NOT [JuJ. 
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The controlled phase-shift gate can itself be considered 
as a derivative of the Ising (zz) coupling gate 
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which corresponds to evolution under the Ising (<J\ z aiz) 
Hamiltonian; Eqns. [2] and ^| are related by single qubit z 
rotations and a global phase shift |T(|. The z rotations 
can be adsorbed into abstract reference frames [|l9| and 
the global phase shift has no physical significance, and 
thus the two equations are essentially equivalent. 

The Ising coupling is of great experimental impor- 
tance, as it provides the basic quantum logic gate for 
many proposed implementations of quantum computing 
24 . Its role is particularly clear in 
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NMR systems, where the scalar coupling (J coupling) 
takes the Ising form in the weak coupling limit |l^] , and 
from this point I will adopt the product operator notation 
widely used in NMR pi 12, However the underly- 



ing principles remain applicable in any system based on 
Ising couplings. 

Consider a system of two spin-i nuclei, I and S. The 
Ising gate (Eq. ||) is implemented by free evolution under 
the J coupling Hamiltonian 



His = 2I z S z 



(4) 



for a time r = (j)/nJ, where J is the coupling strength and 
4> is the desired evolution angle. In order to implement 
accurate controlled phase-shift gates by this means it is 
clearly necessary to know J with corresponding accuracy. 

Obtaining accurate values of J for individual spin-spin 
couplings is easy in small molecules, but is much more 
difficult in larger systems. In particular, consider an ar- 
ray of qubits coupled by Ising interactions, with J cou- 
plings that arc nominally identical but in fact differ from 
one another as a result of imperfections in the lattice. In 
such a system it is desirable to be able to perform some 
accurately known zz evolution over a range of values of 
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J. While this might appear tricky, it is in fact easy to 
achieve using composite pulse techniques. 

The problem of performing accurate zz rotations is 
conceptually similar to that of correcting for pulse length 
errors in single qubit gates fp7| , and the solutions are 
closely related. I begin by describing existing techniques 
for tackling pulse length errors. While a wide variety of 
composite pulse sequences have been developed within 
the context of NMR experiments [H , most of these are 
not suitable for use in quantum computers as they make 
assumptions about the form of the initial state. Instead 
it is necessary to use Class A composite pulses Q which 
work well for any initial state. Two suitable families of 
composite pulses have been described: the SCROFULOUS 
pulse sequences of Cummins et al. [l7| , and the BB1 
family, initially described by Wimperis [p"o]| , and subse- 
quently developed by Cummins et al. [ [l7f . As the perfor- 
mance of the BB1 sequence is much better than that of 
SCROFULOUS I will concentrate on the time symmetrized 
variant of BB1. 

Consider a 9q pulse, which implements a single qubit 
rotation around the x axis on the Bloch sphere jL2|. In 
an ideal system the propagator will be 



U = exp[—i9I x 



(5) 



but in the presence of a fractional error e in the strength 
of the driving field the actual propagator will be 



V = exp[-i6(l + e)I x ]. 



(6) 



The accuracy of this experimental implementation can 
be assessed using the propagator fidelity 



\Tr(yW)\ 
Tr(UW) 
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(note that it is necessary to take the absolute value of 
the numerator as U and V could in principle differ by an 
irrelevant global phase shift). In this case 



T = cos — 
2 
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A better approach is to replace the naive implementation, 
Eq. o, with the composite pulse sequence 
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(9) 



with 4>2 = 30i and cj>i — arccos(— 6 /Air). This gives a 
fidelity expression 

9216 1 ' 

in which the second and fourth order error terms have 
been completely removed JlTj . 

A very similar approach can be used to tackle system- 
atic errors in Ising coupling gates; in effect Ising coupling 



t ^ 4t „ 8t m 4t „ t 
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FIG. 1: Pulse sequence for a robust Ising gate to implement 
a controlled-NOT gate. Boxes correspond to single qubit ro- 
tations with rotation angles of <j> — arccos( — 1/8) ~ 97.2° ap- 
plied along the ±y axes as indicated; time periods correspond 
to free evolution under the Ising coupling, ttJ2I z S z for mul- 
tiples of the time t = 1 /4 J. The naive Ising gate corresponds 
to free evolution for a time 2t. 



corresponds to rotation about the 2I Z S Z axis, and errors 
in J correspond to errors in a rotation angle about this 
axis. These can be parameterised by the fractional error 
in J: 
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Such errors can be overcome by rotating about a sequence 
of axes tilted from 2I Z S Z towards another axis, such as 
2I Z S X . Defining 

6$ = cxp[-i x 6 x (2I Z S Z cos <f) 4- 2I Z S X sin <?!>)] (12) 

allows the naive sequence 9q to be replaced by Eq. || as 
before. The tilted evolutions can be achieved by sand- 
wiching a 2I Z S Z rotation (that is, free evolution under the 
Ising Hamiltonian) between 4> Ty pulses applied to spin S 
|Q. For the case that 8 — tt/2 (which forms the basis 
of the controlled-NOT gate) the final sequence takes the 
form shown in Fig. [l]. 

Clearly it is vital that any robust implementation of a 
quantum gate must be built from components that can 
themselves be implemented robustly. The robust Ising 
gate uses only two components: single qubit rotations 
around the ±y axes, for which robust versions are already 
known jlTj , and periods of evolution under the Ising cou- 
pling. Note that the five time periods must have lengths 
in the integer ratios 1:4:8:4:1, but it is not necessary 
to accurately control the absolute length of the periods, 
as errors in absolute lengths are equivalent to errors in 
the value of J. 

The fidelity gain achieved for coupling gates by this 
approach is identical to that achieved for single qubit 
rotations. In particular for the case that 9 = tt/2 the 
naive fidelity is 



T ^ 1 



32 



while the BB1 approach gives 

^1 ^ 



65536 



(13) 
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Fidelity plots are shown in Fig. |2| Clearly the BB1 ap- 
proach compensates extremely well for small errors in 
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FIG. 2: Fidelity of naive (dashed line) and robust BB1 (solid line) Ising gates as a function of e, the fractional error in J. Plots 
are shown for (a) errors in the range of ±100% and (b) errors in the range of ±10%. Note that in plot (b) the BB1 line is 
indistinguishable from perfection. 



J values, especially within the range ±10%. Over this 
range the infidelity of the BB1 sequence is always less 
that one part in 10 6 . To achieve comparable fidelity with 
a naive gate it is necessary to determine the value of J to 
better than 0.2%. Thus the robust gate can achieve an 
infidelity of 10~ 6 over a range of J values more than 50 
times wider than the naive gate. If even higher fidelities 
are desired the improvement provided by the robust gate 
is even greater. In combination with existing robust sin- 
gle qubit gates, the robust Ising gate provides a complete 
set of robust gates for quantum computation within the 
Ising model. 
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